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Abstract 

Gazeau-Klauder coherent states of the extended trigonometric Scarf potential, underlying the 
quadratic energy spectrum and associated with Jacobi-type X m exceptional orthogonal polyno¬ 
mials Vn ' h '' rn Hx), are constructed. The temporal evolution of wave-packet coherent states are 
performed by means of an autocorrelation function and the full revival properties are investigated 
in the usual time-domain analysis. This latter seems to be less useful for describing the fractional 
revivals due to the complicated nature of coherent wave-packet. Fortunately the autocorrelation 
function revels a little signature of fractional revivals at the vicinity of quarters of the revival time 
T rev due to the quadratic energy spectrum and the use of the wavelet-based time-frequency analysis 
of the autocorrelation function provides an analytical and numerical observation of the fractional 
revivals at different orders of the system. 

PACS numbers: 03.65.-w, 42.50.Ar, 42.50.Md 
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1 Introduction 

Coherent states (CS) were first introduced by Schrodinger [T] in the context of quantum mechanics in 
order to provide a close connection between classical and quantum formalism. They are defined as a 
specific superposition of the eigenstates of the harmonic oscillator in is m and, conventionally, there 
are three ways to define them: 

(HOCS1-) they are a set of eigenstates of the annihilation operator a (Glauber’s approach), 

(HOCS2-) they are regarded as those states obtained by the application of the displacement operator 
D(z) = exp{,W — z*a} upon the vacuum state (Klauder’s approach), and 

(HOCS3-) they minimize the Heisenberg position-momentum uncertainty relation (Schrodinger’s ap¬ 
proach) . 
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Later several proposals for generalizing CS in the context of Lie group structures are suggested 
U El 13 ®, among them we find those constructed by Nieto and Simmons [9], and also by Gazeau 
and Klauder [10]. The former should be considered only in the semiclassical limit which saturate 
a generalized uncertainty relation, while the latter are initially constructed without resort to a Lie 
algebra symmetry and parameterized by two real parameters; i.e. 0 < J < +oo and —oo < 7 < + 00 , 
for models with one degree of freedom with discrete and/or continuous spectra. They are defined by: 


J,l) 


1 y- J n/ 2 exp{-rye n } 

JV(J) vs 


(1.1) 


and p n denotes the moments of a probability distribution p{x) 
nR n 

Pn= x n p(x)dx = TTe;, Po = 1 , ( 1 . 2 ) 

J ° r=i 

where 0 < J < R = lim sup n _ i ._|_ 00 yfp^ and R refers to the radius of convergence. Here e* are the 
spectrum of the Hamiltonian under consideration arranged in a suitable form. 

In Gazeau-Klauder formalism, a coherent state should satisfy the following criteria: (i) continuity 
of labeling, (ii) temporal stability, (iii) resolution of identity, and (iv) action identity. As a consequence, 
several new classes of CS were investigated, such as those often called "nonclassical" which have at¬ 
tracted a lot of attention. Nonclassical effects are characterized by some purely quantum mechanical 
properties m, such as squeezing, photon antibunching, sub-Poissonian photon statistics, and so on. 

However in certain quantum systems with nonlinear (particulary, the quadratic) dependence energy 
spectra, the temporal evolution of wave-packets lead to quantum and fractional revivals na na m. 
Such revivals have been studied extensively in Rydberg atoms [ 15l [16 . T71 iT 8 l 19) . optical paramet¬ 
ric oscillators [201 21], the Jaynes-Cummings model [22, 23], as well as for systems endowed within 
position-dependent mass [231 [23]. Recently, revival dynamics have been put forward, particulary, for 
studying the low dimensional systems [261 EH, the quantum phase transitions [28], the measure of 
nonclassicality by means of Fisher information [29], and also have been related to time-dependent q- 
deformed CS under generalized uncertainty relations m • The full revival arises when the wave-packet 
spreads inside the potential and reconstructs itself during its evolution after a certain time. This phe¬ 
nomenon is observed at long-time, called the revival time T rev (T rev T c j), on which the wave-packet 
is relocalized in the form of a quantum revivals. The fractional revival can be considered as an ad¬ 
ditional temporal structures and found at times equal to rational fractions of the revival time rT rev /l 
with smaller periodicity T c \/l \r/l being an irreducible fraction]. This phenomenon is then subjected 
to partial reconstruction of the initial wave-packet, which are distributed among "l mini-packets" of 
nearly constant value with 1 /l of the total probability density. 

In order to reproduce these revival structures, the autocorrelation function GHEE] is considered 
as widely used technique for describing the revival dynamics of wave-packets and is usually defined by 
the overlap integral 


m 


{^(x,0)\^(x,t)) 

r +00 

/ 'ip*(x : 0)'ilj(x,t)dx, 


(1.3) 


where squared modulus of A(t), |A(t)| 2 , describes the probability of finding the system at time t 
in the vicinity of the initial wave-packet 0). The time-domain analysis is a mathematical tool 
for the autocorrelation function, which use Fourier transform, to reveal perfect quantum revivals for 
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systems whose energy spectrum is purely quadratic in the quantum number |33| . but usually will fail 
to determine precisely the order of the fractional revivals since they appear at specific instants of time. 

To overcome this problem and put into evidence these fractional revivals, an efficient method is 
to apply the continuous wavelet transform (CWT) [341 [35, 36j [33 38j which is powerful for multi¬ 
resolution local spectrum of non-stationary signals. The CWT is a two-dimensional function of the 
scale s and the time shift r referred to as the timescale wavelet representation and equivalent to a 
powerful time-frequency analysis; the scale parameter s stretches (s > 1 ) or compresses (0 < s < 1 ) 
the wavelet till the required resolution while the translation parameter t shifts the wavelet to the 
desired location. The CWT of a function f(t) € L 2 (the vector space of measurable) is defined as 
a decomposition of /(f) into a set of basis functions h SjT (t), often called the wavelets, and generated 
from a simple basic wavelet, known as the mother wavelet, given by 

h s>T (t) = ^=h(*-^-\ (1.4) 

and, using (11.411 . CWT may be regarded as an inner product between the signal, f(t), and a set of 
kernel functions (wavelets) 

W f (s,T) = J f(t)h* jT (t)dt 

= Til mh ’{^ L ) dt (L5) 

However the choice of the kernel function is not specific which represents the main difference 
between the CWT and other well-known transforms such as Fourier, Laplace and Mellin m , and so 
on. Therefore in order to perform calculations, it is important to know what wavelet is used in (11.411 . 

In the present paper, we shall construct coherent states a la Gazeau-Klauder for the Hermitian 
extended Scarf I potential associated with Jacobi-type X m exceptional orthogonal polynomials (EOP). 
These latter are known as the solution of the second-order Strum-Liouville eigenvalues problem with 
rational coefficients [39]. It was seen that although the evolution of the system is effectively governed 
by one time scale, here T rev , its structure at times of fractional revival is not necessarily simple to be 
construct. Thus our main aim of this work is the reconstruction and location of fractional revivals 
corresponding to the Hermitian extended Scarf I potential, by choosing the Morlet wavelet [341 35l [36] 
as a specific kernel function and evaluating the energy density function £(s,t ) = |W|^| 2 (s,r )| 2 for the 
squared modulus of the autocorrelation function, using a wavelet-based time-frequency analysis rather 
than the time-domain analysis of the autocorrelation function. We show that the occurrence and the 
location of fractional revivals for our system is clearly better indicated by CWT than by the standard 
autocorrelation function. 

The organization of this paper is as follows: in Section 2 we describe some generic features of the 
Hermitian extended Scarf I potential in context of X m Jacobi EOP, and construct their associated 
Gazeau-Klauder coherent states. The quantum revival dynamics of the deduced coherent states are 
studied in Section 3 through the construction of their autocorrelation function. In Section 4, we use 
the CWT of the autocorrelation function to reconstruct and localize fractional revivals for the coherent 
wave-packet. Finally, Section 5 is devoted to our conclusion. 
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2 Gazeau-Klauder coherent states for the Hermitian extended Scarf I potential 
associated with X m Jacobi EOP 


In recent work Midya and Roy [39| studied an infinite families of the exactly solvable Hermitian as well 
as non-Hermitian Hamiltonians whose bound-state wave-functions are given in terms of exceptional 
X m Jacobi orthogonal polynomials 


V^’ b ’ m \x) = (-l) r 


Cl + b + j + l _ (_ 0 _i )6 _i) p (a+2,b) ( x 

[ 2(a + j + 1) X ) m [ ’ j ~ l [ ) 


a+ j+l J 


(j = n — m > 0) (2.1) 


with —1 < x < +1, and Vn’ b ’°\x) = Pb a ’ b \x) and P are the standard Jacobi polynomials [40| . The 
real parameters a and b verify the following conditions simultaneously 


(*) b / 0, a, a — b — m T 1 ^ {0,1,... , m — 1} 
(ii) a > m — 2, sgn(a — m + 1) = sgn(6), 

where sgn(a) = a/|a|, or equivalently 


(i) — 1 < b < 0 and m — 2 < a < m — 1 

{ii) b > 0 and a > m — 1. (2.2) 


As a direct application two cases are considered in [39], where one of them concerns the infinite 
numbers of rationally extended real and trigonometric Scarf potential given by: 

C . k 2 (2a 2 + 2b 2 — 1) 2 > k 2 (b 2 — a 2 ) , / 

V ( >(x) = --- -sec kx --- -seckxtankx — 2k m(a — b — m + 1) 

j 2f U , , u , f , , ^ ■ , A p m-l\ sinkx ) 

—k (a — b — m + \ ){a + b + [a — b + \ ) sin kx) —- 

Pm (Sill kx) 


+ 


k 2 (a — b — m + l) 2 cos 2 kx 


Pm-i' > (si n kx) 


(2.3) 


where —ir/(2k) < x < it/(2k) and k ^ 0. Both eigenvalues and eigenfunctions can be obtained 
analytically (with H = 2M = 1) 


£i m) = ^(2n-2m + a + 6 + l) 2 , 


^>(X) = - Sinfc f ^iT!) + Sinfcx)&/2+1/4 V^Hsinkx), 

PX ’ ’(sin kx) 


(2.4) 

(2.5) 


where n = m, m + 1, m + 2 ,... and the normalization constant is given by 




2m + a + b + l)(n — m + a + l) 2 T(n — m + l)T(n — m + a + fe + l) 
2 a + 6+ i(n — 2m + a + l)T(ra — m + a + 2 )T(n + 6 + 1) 


( 2 . 6 ) 


It is to be noted that for m = 0, the potential (12.3|) is reduced to the well-known Scarf I potential. 
The potential given in (12.31) are infinite in number due to the fact that each integer value of m > 0 
gives rise to a new exactly solvable potential and each of them become singular at the zeros of the 
standard Jacobi polynomials, Pm a 1,6 ^(sinfcx), inside the interval — n/(2k) < x < 7r/(2k). These 
singularities can be avoided by applying the restriction (12. 2\ . 

Let us now adapt the material developed above to construct the Gazeau-Klauder coherent states for 
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the extended real Scarf I potential (12.31) . For the case at hand we restrict ourself to the first condition 
of (12.21) . To this end, let us consider the dynamics of the wave-packet in the extended real Scarf I 
potential with eigenvalues ( 12 .4|) : 


, U2 

= — ( 2 n — 2 m + a + b + l ) 2 


= we. 


(m) 


with 


e i m ' ) = ( 2 n — 2 ?n + a + b + iy 


(2.7) 


with the restriction w = k 2 /A. The Gazeau-Klauder coherent states (| 1.1 1) for this system is given by 

+°° rn/2_r -..Am)- 


1 H-oo Tr7 /o r (m)'> 

=<-)(x;J. 7 ) = -±-Y, CXP )l 87en } ^")(x), 

■^( J ) n=0 ^ 


( 2 . 8 ) 


Pn = ( 2 


where 7 = wi. For the coherent states (12. 8 p . p n is easily calculated in terms of Gamma function 
, T(n + 1 + 12 — m) 

r(i + 12 — m) 

where 12 = (a + b + l)/2 and R is infinite. The normalization constant A f(J) is given by 

+ OO 


(2.9) 


jn 

Af’ 2 {J) = — = 1 F ‘2 ( 1 ; 1 + 12 — m, 1 + 12 — m: 


n =0 Pn 


Using the relation m 

r»+oo 


J 


( 2 . 10 ) 


p-r 00 

/ 2J a+p K a _p(2Vj)J n ~ 1 dJ = T(n + 2a)T(n + 2/3), Ren > (-2 Re a, -2 Re /3) (2.11) 

Jo 


we find a solution for the moments of a probability distribution (11.21) 
(J)= A' 0 (\/J) fJ\V-m 


2T 2 (l + n-m) V 4 

where Kq being the modified Bessel function. So, the Gazeau-Klauder coherent states 
become 


s (m) (*;./, 7 ) = 


T(1 + 12 — m) 


J\F 2 ^ 1 ; 1 + 12 — m, 1 + 12 — m; ^ n ~° 


+°° ( J \ n / 2 „— 4 iu(n—m+fl ) 2 t 


T(n + 1 + 12 — m) 


( 2 . 12 ) 

hnally 

(2.13) 


It is easy to see that the Gazeau-Klauder coherent state (12.131) satisfies the following criteria: (i) 
continuity of labeling, (ii) temporal stability, (iii) resolution of identity, and (iv) action identity. For 
example the third condition, associated to the measure AM(J, 7 ), is related to 


r+ir r+oo 

\J,T,m)(J,r,m\dM(J,j) = — J djk( J)| J, 7 ; m)( J, 7 ; m\ dJ, 
if and only if k( J) is defined by 


k{J) = N 2 {J)p{J) = 


k 0 (Vj) 


( J \/ J\ 

—j iF 2 ^1; 1 + 12 - m, 1 + 12 - m; — J 


Jn 


2r 2 (l + 12 — m) 

so the resolution of identity is satisfied. 


(2.14) 


(2.15) 
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3 Full and fractional revivals for extended Scarf I potential 


In this section, we shall review some of the revival dynamics of the coherent states (I2.13|) . It is well 
understood that systems whose energy spectrum is purely quadratic in the quantum number generate 
a perfect revivals [33]■ In our case, although the energy spectra (12.71) are not purely quadratic in n, it 
is a simple matter to redefine the quantum number by shifting n —>• n + m — Q. Then one may suspect 
the existence of perfect quantum revivals for our system. 

For general wave-packet, written in atomic units H = 2 M = 1, 

+oo 

IV’OM)) = ^2c n e~ lEnt \ip(x)), (3.1) 

n =0 


the concept of quantum revivals arises from the weighting probabilities |c n | 2 , with l c n | 2 = 1 - 

So, when the general wave-packets (13.11) play the role of our coherent states (12.131) . then the weighting 
distribution depends on J as 


\Cn{J) | 2 = 


J r 


r 2 (i + n-m) 


j\ n 


(i) 


M 2 {J)Pn 1 F 2 (l;l + fi-m,l + n-m ; ^) r 2 (n + 1 + 12 - m )' 


(3.2) 


In figure 1 we plot the curves |c n (J )| 2 versus n for different values of J. We can see that the 
quasi-Poissonian behavior of weighting distribution is restored at hight values of J and quasi-Gaussian 
shape are almost localized at right gradually as m becomes more and more bigger. In this case, the 




n n 


Figure 1: The weighting distribution given by m for different values of J: (left) a = 2.5, b = 
— 1/2, m = 3 and (right) a = 4.4, b = — 1/3, m = 6 . 


mean values of the number operator N are computed 


(n) 

<n 2 ) 


+oo 

(J, 7 ; m|!V| J, 7 ; m) = ^n|c n (J )| 2 = 

71=0 

+OO 

(J, 7 ;m|lV 2 |J, 7 ;m) =J^n 2 |c n (J )| 2 

71=0 


j \F‘> ^2; 2 + fl — m, 2 + Cl — m; 

4 1 F‘2 + — + U — 

j 2 F 3 ^2,2; 1,2 + U - m, 2 + Q - m\ j'j 
4 iF 2 ^ 1 ; 1 + fl — m, 1 + — m; ^ 


(3.3) 


(3.4) 


in order to display the Mandel parameter Q m defined by 

^ (n 2 ) - (n ) 2 t 

QU = - {n) - L 


(3.5) 
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Figure 2: The Mandel parameter Q m versus J given in (13.51) for different values of parameter a and fixed 
values for b and m: (left) b = — 1/2, m = 4, (center) 6 = —1/4, m = 5 and (right) b = —1/3, m = 6. 


The behavior of CS may be characterized through the Mandel parameter [42] . It is an efficient 
way to characterize non-classical states which have no classical analog. The case of Q m = 0 coincides 
with the definition of coherent states, while the cases of Q m < 0 and Q m > 0 correspond to the 
sub-Poissonian (antibunching effect) and super-Poissonian (bunching effect) statistics, respectively. In 
figure 2 we display the behavior of the Mandel parameter versus J to different states characterizes by 
the parameters a, b, and m. We note that for a given value of a the state starts with super-Poissonian 
behavior, and becomes sub-Poissonian from a certain value of J and thereby maintains this state for 
the rest of the parameter J which hints to a possible squeezing state. We observe, also, that the other 
states start with super-Poissonian behavior (solid curves) but decrease very fast to sub-Poissonian 
(dot-dashed curves) as a increases, then we expect that Q m behaves entirely as a sub-Poissonian state 
from a certain value of a, too. Thus both parameters a and J control entirely the evolution of the 
state. 

As a prerequisite for obtaining quantum and fractional revivals, we suppose that the expansion 
(ED, when the wave-packets are precisely our CS, is strongly peaked around a mean value (n) such as 
the spread An is smaller than (n) ~ n. This allows us to expand the energy eigenvalues (12.71) in terms 
of Taylor series in n around the value n as follows PH 113(111 


fi( m ) ~ fit™'') _L_ 


dE , 


(m) 


dn 


, 1 d 2 E< m > 

(n — n) + 2~^2 


/ —\2 . 1 -PEt" 

{n ~ n) + 6 ^T 


(n - re) 3 H-, (3.6) 


where the derivatives of the energy define distinct time scales, namely, the classical period T c \ = 
2ir/\E^ \', the revival time T rev = 2n/\bE^- l, ' l \", the superrevival time T sr = 2tt/\ bE^ \ r ", and 


so on. 


Neglecting the superrevival time, since the energy ED is a quadratic function in quantum number n 
and considering up to the second-order term, we can write E^ as 


M m) = 4w 


(n — m + Q) 2 + 2(n — m + fl)(n — n) + (n — nY 


which the timescales are define by 
7r 


T c i = 


and 


m 7T 

T — _ 

± rev — „ 

2 ( jJ 


T rev = 2(n — m + fi)T c j. 


(3.7) 


(3.8) 


4 uj (n — m + fi) 

In the present situation, using the last definitions, the Gazeau-Klauder CS reads (in units of 
t = t/T c i) as 


( m )(x t)) = e ~ in ( n -™+ty t 


+oo 


5 >(J)exp{ 


n=0 


—m(n — n)(n + n — 2m + 2f2)_ 
n — m + 


m) 


(*)), (3-9) 
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(a) J=10 


(b) J=20 



Time (In units of T c i) 



Time (In units of T c i) 



Figure 3: The squared modulus of the autocorrelation function plotted from (13.1211 along t = 2 T rev for 
a = 4.4, b = —1/3 and m = 6. The time t is given in units of T c \. 


where c n (J ) are defined by (13.2(1 . For the coherent states (13.911 it is possible to calculate the quantum 
revivals through the autocorrelation function defined as the projection of a wave-packet CS onto its 
initial state 


A(t) = <S (m >(z,0)| 




i F 2 (1; 1 + n - m, 1 + 12 — ?n; 4 ) n = 


=o 


r 2 (?r — m + 12 + 1) 


Using the relation 


+oo 

E' 


+oo 


T. anxn ) ( ^ 


+oo 


Cn — 


^ ^ & kbn—k ? 


with 

n =0 v n=0 7 v n=0 7 k=0 

the straightforward calculations lead to express the squared modulus of A(t ) as 


\m \ 2 = 


r 4 (i + n - m) 


+oo n 

EE 


exp 


f «7r(n— 2k)(n—2m+2Q)j\ ( J\i 
\ n-m+Q l 4/ 


(3.10) 


(3.11) 


i Ff (l; 1 + Q - m, 1 + n - m; (\ ^ m + + 1 ) I ’ 2 ( n - ^- ^ + 0 + 1) 


( 3 - 12 ) 


where in order to make numerical calculations we need to truncate the series at some order in n, 
keeping in mind that n > m. 

The evolution of \A(t = f/T c i)| 2 is shown in figure 3 in the cases J = 10, 20, 40 and 100. These 
plots are generated by choosing n max = 50 and n = 100 for a = 4.4, b = —1/3 and m = 6. The 
revival time calculated through (13.121) is T rev — 2896.825 (in units of T c \) and using (13.81) one find 
T c i ~ 15. All cases clearly show that the time evolution of the system is periodic with period T rev , 
which is an exact revival time too, since the revival structures are not modulated by a superrevival 
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time. The sharp peaks appearing in the cases J = 20, 40 and 100 are characterize by the occurrence of 
fractional revivals which become more and more apparent as J increases. The function (I3.12p proves to 
be symmetric about T rev /2 which means that all the features of coherent wave-packet evolution appear 
at a time interval from 0 to T rev /2. After that time all the events will repeat again and again. As we 
can see in figures 3(c) and (d) especially interesting are fractional revivals occurring at quarters of the 
revival time (indicated by an arrow) showing a suspect symmetry which indicates a strong occurrence 
of fractional revivals at that times since the energy is quadratic in n. Unfortunately these peaks 
are uniformly mixed so that we do not have any idea on which value of the autocorrelation function 
appears at what time; thus fractional revivals are still invisible and should be resolved. Then a natural 
question arises: what about the reconstruction and location of fractional revivals? This question will 
be considered in more details in a forthcoming section. 


4 Reconstruction and location of fractional revivals from wavelet-based time- 
frequency representation 


In this section our key aim is the observation of fractional revivals of the potential m since the 
autocorrelation function fails to reconstruct them. An efficient and especially illustrative method 
for this kind of construction involves the use of the time-frequency analysis based on the continuous 
wavelet transform (CWT) [54} [35, 36]; this latter is basically a two-dimensional function of the scale 
s (inversely to the frequencies) and the time shift r. The specific choice of using CWT rather than 
short-time Fourier transform (STFT) consists in the fact that STFT uses a single analysis window, 
while CWT uses short windows at high frequencies and long windows at low frequencies. This property 
overcomes the limitation of STFT and allows us to uncover finer details of signal acting at smaller 
scales. 

Here, in our case, the signal is our autocorrelation function given by (13.121) and the kernel function 
is the well-known Morlet wavelet defined by 

h(t) = -^ JI e i ^ t e- t2 / 2 , (4.1) 

where ujq is called the central frequency. Let us now be more specific by redefining t —>• £ = (t — t)/s, 
the CWT of |A(f)| 2 given in (11.51) can be written as 

/ +oo 

\A(s( + T)\ 2 h*(()d(, (4.2) 

-OO 

where h(Q = 7r~ 1 / 4 e^ 0< ’e _ ^ 2 / 2 and by performing the integration over £ we get 


W\a\2(s,t) = 7r 1 / 4 v / 27 ts EE \c n (J)\ 2 \ck(J)\ 2 exp | - 2m(n- ^—— Ir 

( l\ J 1 n + k — 2 n\l 2 i , . 

X eXP 1 2 W0 + 27TS ^ U ~ k \T~i + - T~ -) ( 4 - 3 ) 


Despite we do not know an exact expression for the central frequency ujq, it is useful to identify, at 
first approximation, ojq ~ u = k 2 / 4 in order to calculate the minima of LF|^| 2 (s,t). For convenience, 
we reproduce here, in part, some of results obtained by Ghosh and Banerji m- Since the weighting 
functions |c n (J)| 2 and |cfc(J)| 2 appearing in (14.31) are peaked around ri. the central frequency can be 
calculated by replacing n by n and k by n + p, where p is a positive integer. This gives rise to the 
relation 


uj = 


2nsp f 


Tel 


+ P - 


T r 


cl 


27TS 


P, 


Tel 



(4.4) 
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and since the scale parameter s is related to frequency v by the relation s = ui/(2m'), then we obtain 
the formulae 


v ~ 


P 


and 


u 


s = 


T c i Ittu 

Each frequency corresponds a particular time r given by 


(4.5) 



(4.6) 


where q is a positive integer. If the restriction q < 2p holds, then one can see the correspondence 
between (14.611 and the well-defined fractional revivals of the autocorrelation function Aft) localized at 
various rational multiples of the revival time given by t = ( r/l)T iev (mod l). Both expressions (14.51) 
and (14.6[) provide the required resolution in frequency and the desired location in time for the recon¬ 
struction and location of fractional revivals; in other words, one can find both time and frequency in 
terms of the two timescales T rev and T C 1. 

Ghosh and Banerji used their results to correctly analyze the localization of different patch cor¬ 
responding to the fractional revival it contributes to, using the Time-Frequency ToolBox for MatLab 
[35] , Here we test the robustness and efficiency of CWT as a mathematical filter, in order to filter out 
fractional revivals from the complicated plots of figure 3 obtained from the autocorrelation function. 
To this end we restrict our analysis to the reconstruction and localization of fractional revivals. We 
illustrate this by plotting in figure 4 the squared modulus of both Aft ) and W\j^i(t). The figure (4a) 
displays the evolution of |A(i)| 2 given by ()3.12D in units of T c \ along t = 10T rev , while figure (4b) gives 
the time dependence of |IH|yi| 2 (r)| 2 given by (14. 3 p in units of T rev for the first harmonic fp = 1) and at 
a time interval from 2 T rev /9 to 2 T rev /7, i.e., near the first quarter revival t — T vey /4:. As shown on the 
figure (4b), the coherent wave-packets evolve in time and split into a collection of spatially distributed 
mini-" empire state building" groups, each of them closely reproduces the shape of the initial coherent 
wave-packets of figure (4a), and are all of nearly constant values. The motion of these fractional re¬ 
vivals is periodic with period given by a rational fractions of T c \ ~ 15 calculated from (13.811 . as can 
be seen on the figure (4b). Thus, for our system under study, we can conclude that the wave packet 
remains coherent for a number of classical periods. 

Now to correctly analyze the location of fractional revivals, we use m and m to express 
VT |^|2 (s,r) in terms of the parameters p and q, which gives 


W\ A \2(p,q) = 


/ rji “bOO H“00 z' / fj-i 

'-"■/-EE |c n (J)| 2 |c fc (J)| 2 exp i - 2mfn - k) 


p 


x exp < — 


n= 0 k =0 v 

u 2 T n — k ( , T c i 

1 H-( 1 + (n + k - 2 n) —— 

P V Hey 


m + n + k — 2n ) — 
Tel ) 2 P 


(4.7) 


and we display in figure 5 the CWT of the autocorrelation function (14.3p and the density plot of 
lndWi^p^, (/)| 2 ), for different values of p and q. The figure (5a) shows how CWT of the autocorrelation 
function Aft) displays a possible location of fractional revivals at both absolute minimum (vertical 
arrow) and relative minima (dashed arrows) of W^^fr), for the first four harmonics p = 1, 2,3, and 4. 
The absolute minimum clearly predicts that the first quarterly fractional revival should occur around 
T = r/T rev ~ 1/4, and fade away for r > 1/4. One can observe that they are more better localized 
than those occurring at vicinity of r ~ 3/4. These imperfect and asymmetric quarter revivals take the 
slightly modified form due essentially to the fact that the energy spectrum are not purely quadratic 
in the quantum number n. However, the relatively minimum areas can be the scene of several other 
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(4b) 
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Figure 4: (4a) The squared modulus of the autocorrelation function A(t) defined by (13.1211 along 
t = 10T rev , (4b) the squared modulus of W|A| 2 ( r ) defined by (14.3() for p = 1. Both figures are plotted 
for parameters: a = 4.4, b = —1/3, m = 6 and J = 10. 

orders of fractional revivals but with a less significant probability amplitude than those detected in the 
absolute minimum areas. 

Indeed, referring to the figure (5b), one can observe that the occurrence of fractional revivals 
is clearly better indicated by ln(|VFMi 2 (p, q)\ 2 ) than by |A(t)| 2 . A characteristic straight dark rays 
are clearly visible on the p — q plan, which correspond to the different areas (patches) where the 
fractional revivals occur and denoted by the different colored points. Each point is well-localized about 
a particular p (i.e., frequencies v) and q (i.e., time r) and located at the center of the corresponding 
patch; it tells us which the fractional revival it contributes to. For examples, using (14.61) . we can observe 
that the red points are all localized at the first quarter revival, i.e., t = T rev /4, while the corresponding 
value of time associated to the yellow points is localized at the second one, i.e., at t = 3T rev /4. Thus 
the evolution of the system shows clearly the occurrence of fractional revivals exactly at quarters of the 
revival time, with a little imperfect and asymmetric modified form of rays due to the raison mentioned 
above. However the pink and purple points, both localized on the white regions, are all centered at 
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( 5 a) 


(5b) 



Figure 5: (5a) The plot of the CWT of the autocorrelation function defined by (14.31) in units of 
r = r/T rev , (5b) log density plot of the squared modulus of (14.71) . 


times r = T rev /2 and r = T rev , respectively. This indicates clearly that one does not observe the 
occurrence of the fractional revivals at these times. 

On the other hand, the existence of relatively minimum areas clearly manifest a slight emergence 
of a few patches showing the location of fractional revivals at the time r = T rev /8 (black points), 
r = T rev /6 (blue points), and r = T rev /3 (green points), respectively, and disappeared completely in 
the interval T rev / 2 and T rev , except of course for those associated to t = 3T rev /2 (yellow points). These 
occurrences are restrictive in the sense that are localized in relatively minimum areas, as displays in 
figure 5(a), with a certain probability densities less than those attributed to patches located at exactly 
quarters of the revival time. 

5 Conclusion 

In this paper, we have constructed the Gazeau-Klauder CS for the extended Scarf I potential associated 
with exceptional X m Jacobi OP. These CS satisfy the criteria of continuity of labeling, resolution of 
unity, temporal stability, and action identity. The probability distribution shows a Gaussian behavior 
for a large values of J. The Mandel parameter starts with super-Poissonian behavior and decreases very 
fast to sub-Poissonian as a and J increase and remains in this state for the rest of these parameters, 
which indicates that our CS (I3.9[) exhibit squeezing. 

Using the autocorrelation function of (13.911 . a full revivals are obviously occur at each multiple 
of T rev and sharp peaks revel a little signature of the fractional revival structure. These latter are 
uniformly mixed, so that their orders can not be determined exactly. To overcome this problem, 
we have made used of the continuous wavelet transform of A(t) in order to reconstruct and localize 
fractional revivals. This study has shown (i) that fractional revivals faithfully reproduce the shape 
of the initial coherent wave-packets, within an imperfect and asymmetric quarterly form due to the 
fact that the energy spectrum are not purely quadratic in n, and (ii) the location of fractional revivals 
is clearly better viewed by ln(|kFMi 2 (p, q)\ 2 ) than by |yl(f)| 2 , which indicates that the wavelet-based 
time-frequency representation can be regarded as a powerful tool and a complementary method for 
reconstructing and localizing fractional revivals. 
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